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We present a theoretical study on the response properties to an external electric field of strongly 
correlated one-dimensional metals. Our investigation is based on the recently developed Bethe- 
Ansatz local density approximation (BALDA) to the density functional theory formulation of the 
Hubbard model. This is capable of describing both Luttinger liquid and Mott-insulator correlations. 
The BALDA calculated values for the static linear polarizability are compared with those obtained 
by numerically accurate methods, such as exact (Lanczos) diagonalization and the density matrix 
renormalization group, over a broad range of parameters. In general BALDA linear polarizabilitics 
are in good agreement with the exact results. The response of the exact exchange and correlation 
potential is found to point in the same direction of the perturbing potential. This is well reproduced 
by the BALDA approach, although the fine details depend on the specific parameterization for the 
local approximation. Finally we provide a numerical proof for the non-locality of the exact exchange 
and correlation functional. 

PACS numbers: 



I. INTRODUCTION 

Material systems, whose electronic structure cannot be 
described at a mean field level, are conventionally named 
strongly correlated. These display an enormous variety 
of properties, which all originate from the interplay be- 
tween Coulomb repulsion and kinetic energy, and from 
their dimensionality. Phenomena related to electron- 
electron correlation include metal-insulator transition, 
Tomonaga-Luttinger liquid behaviour and superconduc- 
tivity, just to name a few [IJH] . In particular electron cor- 
relations play a fundamental role in one-dimension (ID). 
In ID confined structures, electrons cannot avoid each 
other and collective excitations emerge over the ground 
state, so that the Fermi liquid picture breaks down. In 
fact one can demonstrate that the ground state of an in- 
teracting ID object is always a Luttinger liquid regardless 
of the strength of the electron-electron interaction [TJ [5] . 
Although some aspects are still controversial, experimen- 
tal evidence supporting the existence of Luttinger liquids 
in ID has been provided for carbon nanotubes [3] and for 
atomic wires built of surface terraces [H [5] . 

Strongly correlated systems are regularly modeled by 
means of effective Hamiltonians, which usually lack all 
the details of an ah initio description, but capture the 
relevant physical properties arising from electron correla- 
tion. The advantage of dealing with effective Hamiltoni- 
ans is that they are commonly mathematically tractable 
and general enough to be applied to a variety of prob- 
lems. Among the many effective Hamiltonians that one 
can construct the Hubbard model [SHE] has enjoyed a 
vast popularity since it is simple and still can capture 
the subtle interplay between Coulomb repulsion and ki- 
netic energy. 

Although exact solutions of the Hubbard model are 
known in particular limits [9], a general one for an ar- 
bitrary system, which can be finite and inhomogeneous, 



requires a numerical treatment. This however represents 
a severely demanding task, since the Hilbert space asso- 
ciated to the Hubbard Hamiltonian for L sites is 4 L di- 
mensional, so that exact (Lanczos) diagonalization (ED) 
can only handle a relatively small number of sites. Other 
many body approaches, such as the density matrix renor- 
malization group (DMRG) [lOl [11] , extend the range to 
a few hundred sites, but little is possible beyond that 
limit. It would be then useful to have a method capable 
of describing accurately the ground state and still having 
the computational overheads of a mean field approach. 
Such a method is provided by lattice density functional 
theory (LDFT). 

LDFT was initially proposed by Gunnarsson and 
Schonhammer [321 03] as an extension of standard, ah 
initio, DFT P31 |T5] to lattice models. The theory es- 
sentially reformulates the Hohenberg-Kohn theorem and 
the Kohn-Sham construction in terms of the site occupa- 
tion instead of the electron density. Although originally 
introduced with a pedagogical purpose, LDFT has en- 
joyed a growing success and it has been already applied 
to a diverse range of problems. These include funda- 
mental aspects of DFT and of the Hubbard model, as 
the band-gap problem in semiconductors Q2] , the dimer- 
ization of ID Hubbard chains [TB] and the formation of 
the Mott-Hubbard gap 17J. LDFT has also been em- 
ployed for investigating effects at the nanoscale traceable 
to strong correlation, like the behavior of impurities [TS] , 
spin-density waves [115] and inhomogeneity [50], as well 
as more exotic aspects like the phase diagram of har- 
monically confined ID fermions [2l] and that of ultra- 
cold fermions trapped in optical lattices [22H23]. More 
recently LDFT has been extended to the time-dependent 
domain [25], to quantum transport [26] and to response 
theory [57]. 

As in standard DFT also LDFT is in principle ex- 
act. However its practical implementation is limited by 



2 



the accuracy of the unknown exchange correlation (XC) 
functional, which introduces the many-body effects into 
the theory. The construction of an XC functional begins 
with choosing a reference system, for which some exact 
results are known. These impose a number of constraints 
that the XC functional must satisfy, as for example its 
asymptotic behavior or its scaling properties. Then the 
functional is built by interpolating and fitting to known 
many-body reference results. Such a construction for in- 
stance has been employed in the case of the local density 
approximation (LDA) in ab initio DFT. The reference 
system in two and three dimensions is usually an elec- 
tron gas of some kind, since one aims at reproducing a 
Fermi liquid. However in ID the known ground state has 
a Luttinger-liquid nature, so that the reference system 
should be chosen accordingly. In the case of the Hubbard 
Hamiltonian in ID a powerful result is that obtained by 
Lieb and Wu [28] for the homogeneous case by using the 
Bethe Ansatz. This is the basis for constructing an XC 
functional for the Hubbard model in ID [TT1 118j . 

In this work we evaluate the ability of a range of known 
approximations to the XC functional for the ID Hubbard 
model at predicting the electrical response to an external 
electric field of finite ID chains away and in the vicinity 
of the Mott transition. This is relevant not just as a test 
for Hubbard LDFT but also for understanding real mate- 
rials, whose electrical response can be mimicked in terms 
of the Hubbard model [351130]. The ID case in particular 
can provide important insights into the nonlinear optical 
properties of polymers [31) . Our strategy is that of con- 
stantly comparing the DFT results with those obtained 
with highly accurate many-body schemes. For these we 
use exact diagonalization for small chains and the DMRG 
method for larger systems. Our calculations reveal a sub- 
stantial good agreement between LDFT and exact results 
for both the polarizability and the XC potential response 
of finite ID chains. The paper is organized as follows. 
In the next section we will briefly review the Hubbard 
LDFT and the approximations used for constructing the 
XC functional. Then we will discuss results, first for the 
electrical polarizabilitics and then for the response of the 
XC potential to an external electric field. Finally we will 
carry on a numerical investigation on the validity of the 
local approximation to the XC functional and then we 
will conclude. 



II. THEORETICAL FORMULATION 

One-dimensional correlated metals can be described by 
the homogeneous Hubbard Hamiltonian, H\j. For a ID 
chain comprising L sites H\j writes 

L-l L 

Hy = -t ^2 (c\ +la c ia + he) + *Wtra»4-> (!) 

1=1, a i—1 

where the first kinetic term describes the hopping of elec- 
trons with spin a {a = j" , I) between nearest neighbour 



sites with amplitude t > 0, while the second accounts for 
the electrostatic repulsion U > of doubly occupied sites. 
In the equation |l| c\ a {ci a ) is the fermion creation (anni- 
hilation) operator for an electron at site i with spin er and 
the site occupation operator is written as h ia = c\ a Ci a . 
Clearly there is only one energy scale in the problem so 
that the ratio U/t determines all the electronic proper- 
ties. Note that a second energy scale can be included 
in the problem by adding to the Hamiltonian an on-site 
energy term Y]f—i „ Cih ia mimicking a ionic lattice. 

As discussed in the introduction the fundamental 
quantity of LDFT is the site occupation, m , which is cal- 
culated by solving the equivalent Kohn-Sham problem. 
This can be generally written as 

L 

i + h-xi) + vk s ]^ a) = 6 (Q M Q) , (2) 

where v^ s is the general Kohn-Sham potential. The oc- 
cupied Kohn-Sham eigenvectors, 4>[ a \ define rij 

* = £™ (a) i*, (a) i a , ( 3 ) 

a 

where are the occupation numbers, which satisfy 

J] w( a ' = N with N being the total number of elec- 
trons. By following in the footsteps of standard ab initio 
DFT the Kohn-Sham potential can be written as the sum 
of three terms 

"ks = K + *4t + «xc] . ( 4 ) 

where «h = Urii/2 is the Hartree potential and i>g Xt is 
the external one. The last term in equation Q is the 
XC potential, which needs to be approximated. 

The Kohn-Sham equations simply follow by variational 
principle from the minimization of the energy functional. 
Thus the total energy of the system, E, can be defined 
as 

E[{ ni }} =J2 w{a) * {a) -T, v xcn-Yl ~+£xc[K}], 

a i i 

(5) 

where the last term is the XC energy. Note that different 
values of U and t define completely the theory, so that 
one has a different functional for every value of U/t. 

We now review the strategy used for constructing a 
suitable local v^ c [17U19] . The guiding idea is that of 
defining v^ c as the local counterpart of the Bethe Ansatz 
potential for the homogeneous Hubbard model (for an 
infinite number of sites), i.e. 

i&c|balda = «xc m (M, U)\ n ^ ni . (6) 

Here BALDA stands for Bethe Ansatz local density ap- 
proximation and <° c m (n, t, U) is the XC potential for the 
homogeneous Hubbard model, which is defined only in 
terms of the band filling n = N/L (note that for the ho- 
mogeneous case n, = n for every site i). Formally, and 
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in complete analogy with ab initio DFT, v^^(n,t, U) is 
obtained by functional derivative of the exact energy den- 
sity, e(n,t,U), of the reference system (in this case the 
homogeneous Hubbard model), after having subtracted 
the kinetic energy density of the non-interacting case 
e(n,t,U — 0) and the Hartree energy density, en(n, U). 
This gives us 

d 



v^(n,t,U) = ^[e(n,t,U)-e(n,t,U = 0)-e u (n,U)} . 

(7) 

The question is now how to obtain e(n,t, U). Two al- 
ternative constructions have been proposed in the past 
and here we have adopted and numerically implemented 
both. The first one consists in using the analytical pa- 
rameterization proposed by Lima et al. [TH] [19] , which 
interpolates the known exact results for: 1) U — > and 
any n < 1, 2) U — > oo and any n < 1 and 3) n = 1 and 
any U. The resulting XC potential can then be written 



hom 



(n, t, U) = tfX 
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, (8) 



where k = 1 — \n — 1|, /J, = sgn(n — 1) and (3(U) is a in- 
dependent parameter, which can be determined by solv- 
ing a transcendental equation. The alternative route is 
that of employing a direct numerical solution of the cou- 
pled Bcthe Ansatz integral equations. This approach has 
been already used for the study of ultracold repulsive 
fermions in ID optical lattices [2"2"] . 

The first parameterization is known as BALDA/LSOC 
[3"2] and the second as BALDA/FN (FN = fully numer- 
ical) . In figure [TJ the XC potentials as a function of the 
electron filling for the two schemes are shown for differ- 
ent values of U. In the picture (and in the calculations) 
we always use the particle-hole symmetry, which imposes 

^xc m ( n > M» U) = -w X c m ( 2 - n > *> u )- From the n S ure 
one can immediately observe that the potential in both 
cases has a discontinuity in the derivative at half-filling 
(n = 1, N = L). This reflects the fact that the underlying 
homogeneous ID Hubbard model has a metal-insulation 
transition for n = 1. Such a discontinuity in the deriva- 
tive of the potential, as in standard ab initio DFT, is re- 
sponsible for the opening of the energy gap. The second 
observation is that the two parameterizations always co- 
incide by construction at n = and n = 2 but that their 
agreement over the entire n range depends on the value 
of U. In particular one can report a progressively good 
agreement as U increases. This is not a surprise since 
the BALDA/LSOC potential is constructed to exactly 
reproduce the U — > oo limit. 



III. POLARIZABILITIES 

We calculate the electrical polarizability of linear 
chains with the finite difference method, i.e. as numerical 
derivative of calculations performed at different external 




FIG. 1: (Color online) Exchange-correlation potential 
VxcXM,^) of the ID Hubbard model as a function of 
the electron filling, n, for different values of interaction 
strength, U. Here we report data for both BALDA/LSOC 
and BALDA/FN. Note that the agreement between the two 
schemes improves as U increases. 



electric fields. An external electric field enters into the 
problem by adding to the Hubbard Hamiltonian H\j the 
term 



e£x — e£ — x)c\ci , 



(9) 



where x — g(Z> + 1) is the middle site position of the 
chain, e is the electronic charge (e = —1) and £ is the 
electric field intensity (the electric field is applied along 
the chain). In general the electrical dipole, P, induced by 
an external electric field can be calculated simply as the 
expectation value of the dipole operator over the ground 
state wave-function |^o(^)) (note that this is a general 
definition so that \^ (£)) is not necessarily the Kohn- 
Sham ground-state wave-function), i.e. 



P = e{$ (£)\J2(i-x)4c i \$ (£)) 



dE (£) 
d£ ' 



(10) 



where Eq is the ground state energy. For small fields P 
can be Taylor expanded about £ = so that the linear 
polarizability, a, is defined as 



P - a£ + - ( £ 3 + 0(£ 5 ) 



d 2 E Q (£) 
d£ 2 



(11) 



Our calculation then simply proceeds with evaluating 
E (£) for different values of £ and then by fitting the 
first derivatives with respect to the field to the equa- 
tion |TT] ), as indicated in reference [31]. We note that 
our finite difference scheme is not accurate enough for 
calculating the hyper-polarizability, 7, which then is not 
investigated here. 
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It has been already extensively reported that BALDA- 
LDFT gives a substantial good agreement with exact cal- 
culations in terms of ground state total energy [T7J [TB] . 
The polarizability however offers a more stringent test for 
the theory since it involves derivative of Eq . Hence it is 
important to compare the various approximations with 
exact results. For small chains, L < 18, these are ob- 
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and consequently a reduction of a. Away from U 
the dependence of a on U ft can be fitted with 
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FIG. 2: (Color online) Linear polarizability, a, as a func- 
tion of the Coulomb repulsion U/t. Results are presented 
for BALDA/LSOC and BALDA/FN and they are compared 
with those obtained with either exact diagonalization (ED) 
or DMRG calculations. In the various panels we show: (a) 
L — 12 at quarter filling (n = 1/2), (b) L = 16 at quarter 
filling, (c) L = 60 and N = 20, and (d) L = 60 at quarter 
filling. 

tained by simply performing ED. However for the longer 
chains ED is no longer feasible and we employ instead 
the DMRG scheme [HI US]. DMRG has been widely used 
to investigate one-dimensional and quasi one-dimensional 
quantum systems. It usually performs best with open 
boundary conditions and utilizes appreciable computa- 
tional resources depending on the number of states that 
are kept for the calculation. Our DMRG calculations are 
performed by employing the Algorithms and Libraries 
for Physics Simulations (ALPS) [31] package for strongly 
correlated quantum mechanical systems. The DMRG re- 
sults are obtained by using a cutoff of m = 350, i.e. by 
retaining the dominant 350 density matrix eigenvectors. 

Let us start our analysis by looking at the polarizabil- 
ity as a function of the energy scale U/t. Selected results 
for quarter-filling, n = 1/2, and for n — 1/3 are presented 
in the various panels of figure [2] Note that throughout 
this work we always stay away from the half-filling case 
(n = 1), where the derivative discontinuity of the poten- 
tial makes the LDFT convergence problematic. 

In general we find that the polarizability decreases 
monotonically with increasing the on-site repulsion U. 
This is indeed an expected result since an increase in on- 
site repulsion means a suppression of charge fluctuations 



a(U /t; L, n) = ao(L, n) 



(12) 



where all the parameters have a dependance on the length 
of the chain and on the band filling. The results of such 
a fitting procedure are reported in table [I] Note that in 
the fit we did not impose any constraints and we have 
included only points with U/t> 1. 



Method 


L N n ao £ 


ED 


12 6 1/2 59.69 0.23 


balda/lsoc 


62.06 0.27 


balda/fn 


59.50 0.25 


ED 


16 8 1/2 142.88 0.27 


balda/lsoc 


135.07 0.31 


balda/fn 


143.56 0.30 


DMRG 


60 30 1/2 8939.5 0.32 


balda/lsoc 


8673.1 0.33 


balda/fn 


8837.8 0.31 


DMRG 


60 20 1/3 6931.6 0.29 


balda/lsoc 


6401.0 0.30 


balda/fn 


6920.7 0.29 



TABLE I: Scaling parameters for a(U/t;L) as obtaine d by 
fitting the data of Fig. [2] to the expression of equation ( 12 1 . 
Note that the fit has been obtained without any constraints 
and by including data only for U/t> 1. 



From the fit and from figure [2] one can immediately 
note that both the BALDA flavors of the exchange and 
correlation functional reproduce rather well the exact re- 
sults, in good agreement with previously published cal- 
culations [27] . The agreement is particularly good for the 
FN functional, which matches the ED/DMRG results al- 
most perfectly over the entire range of U /t's and filling 
investigated. A quantitative assessment of goodness of 
the BALDA results is provided in figure [3] where the rel- 
ative error, 8, from the reference exact calculations is 
presented. In general, and as expected, we find that the 
error grows with U/t, i.e. with the system departing from 
the non-interacting case. However, there is also a satu- 
ration of the error as the interaction strength increases, 
reflecting the fact that both the BALDA potential are 
exact in the limit of U — > oo. As a further consequence 
of the U — > oo limit, we also observe that the relative 
error between BALDA/LSOC and BALDA/FN reduces 
as U grows. 

Given the accuracy of the BALDA/FN scheme we have 
decided to use the same to investigate in more details 
the scaling properties of a(U/t;L). First we look at the 
scaling as a function of the interaction strength U/t. In 
this case we always consider a chain containing L = 60 
sites for which the deviation from the DMRG results is 
never larger than 2%. Furthermore this is a length which 
allows us to explore a rather large range of electron filling, 
so that it allows us to gain a complete understanding 
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FIG. 3: (Color online) Relative error between BALD A calcu- 
lated polarizabilities and those obtained with exact methods 
(either ED or DMRG). In the panels we show: (a) L — 12 
at quarter filling (n = 1/2), (b) L — 16 at quarter filling, (c) 
L = 60 and N = 20, and (d) L = 60 at quarter filling. 




FIG. 4: (Color online) Polarizability as a function of U/t 
for a chain of 60 sites and various filling factors, n. The 
figure lege nd reports the fitted values for the exponent £ [see 
equation ( 12 l] . The symbols represents the calculated data 



while the solid lines are just to guide the eyes. In the inset we 
present the exponent £ as a function of the filling factor n. 



of the scaling properties. Our results are presented in 
figure [4] where we show a as a function of U jt for different 
filling factors, we list the values of £ obtained by fitting 
the actual data for U/t > 1 to the expression in equation 
(12) and we provide (inset) the dependence of £ on n. 

In general the fit to our data is excellent, suggesting 
the validity of the exponential scaling of the polarizability 
with the interaction strength (away from half filling). In 
particular we find that £ decreases monotonically with n 
for n > 0.2 but it increases for smaller values. This means 



that £(n) has a maximum just before n = 0.2, which 
appears rather sharp (see inset of figure [4]). We are at 
present uncertain about the precise origin of such a non- 
monotonic behavior. However, as we will see in details 
later on, we notice that the response of the exchange 
and correlation potential to the external electric field has 
an anomaly for small U and n. We believe that such 
an anomaly might be the cause of the non-monotonic 
behaviour of £. 

Next we turn our attention to the scaling of a with the 
chain length. In figure[5]we present a(L) for two different 
filling factors (n = 1/3 and 1/2) and different values of 
U/t. Data are plotted both in linear and logarithmic 
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FIG. 5: (Color online) Scaling of the polarizability as a func- 
tion of the chain length, L. Panel (a) and (b) are for n — 1/3 
while (c) and (d) for n — 1/2. Note the linear dependence of 
the a(L) curve when plotted on a log-log scale, proving the 
relation a(L) = aiL 1 

scale, from which a clear power-law dependance of a on L 
emerges. A fit to our data provides the following scaling 



a(U/t;L) = a^V 1 



(13) 



Importantly this time we find essentially no dependance 
of both a.\ and 7 on either U jt or n. The fit reveals 
a value for the exponent of 7 ~ 3 (the range is from 
7 = 2.93 to 7 = 2.98). This is what expected for free 
electrons in ID [31] , and it is substantially different from 
the predicted linear scaling at n = 1. Our results thus 
confirms that away from n = 1 the electrostatic response 
of the Hubbard model is similar to that of the non- 
interacting electron gas. Going in more details we find 
a rather small monotonic dependance of 7 on U jt. This 
however depends also on n since for n — 1/3 we find that 
7 reduces as U/t is increased (from 2.98 for U/t = 0.5 
to 2.93 for U/t = 100), while the opposite behavior is 
found for n = 1/2 (7 = 2.94 for U/t = 0.5 and 2.96 for 
U/t = 100). 
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IV. RESPONSE OF THE BALDA POTENTIAL 
TO THE EXTERNAL FIELD 

In ab initio DFT the failures of local and semi-local 
XC functionals in reproducing accurate linear polariz- 
abilities are related to the incorrect response of the XC 
potential to the external electric field |35j [36] , which in 
turn originates from the presence of the self-interaction 
error (371 I38j . In particular for ab initio DFT the ex- 
act XC potential should be opposite to the external one, 
while the LDA/GGA (generalized gradient approxima- 
tion, GGA) returns a potential which responds in the 
same direction. In order to investigate the same feature 
for the case of the Hubbard model LDFT we calculate 
the potential response 

where «xc( n i) * s the exchange and correlation potential 
at site i in the presence of an electric field £ . Also in this 
case we adopt the finite difference method and we use 
£ = 0.01, after having checked that the trends remaining 
unchanged irrespectively of the field strength. 

In order to provide a benchmark for our calculations 
we also need to evaluate the potential response for the 
exact Hubbard model. We construct the exact poten- 
tial by reverse engineering, a strategy introduced first by 
Almbladh and Pedroza [35] and by von Barth [ID] and 
then applied to both static and time dependent LDFT 
by Verdozzi [25] . This consists in minimizing about the 
Kohn-Sham potential the functional J- (in reality here 
this is just a function) defined as 

L 

J>xc]=l>f s -nr act ) 2 , (15) 

i 

where ri° xact is the exact site occupation at site i as ob- 
tained by either ED or the DMRG method, while nf s is 
the Kohn-Sham one. 

Our results are summarized in figures [6] and [7] where 
we show A«xc as a function of the site index for a 60 
site chain occupied respectively with 10 (n = 1/6) and 
30 (n = 1/2) electrons. The external electrostatic poten- 
tial here decreases as the site number increases, i.e. it 
has a negative slope. Results are presented for DMRG, 
BALDA/LSOC and BALDA/FN and for different values 
of U/t. 

In general and in contrast with ab initio DFT, we find 
that the response of the exact Hubbard-LDFT XC po- 
tential is in the same direction of the external pertur- 
bation for both the filling factors investigated and re- 
gardless of the magnitude of U/t. The response however 
becomes larger as U/t is increased (the slope of Avxc 
is more pronounced), a direct consequence of the fact 
that for large U's small deviations from an homogeneous 
charge distribution produce large fluctuations in the po- 
tential. Such a behaviour is well reproduced by both the 
BALDA functionals, with the BALDA/FN scheme per- 
forming marginally better than the BALDA/LSOC one, 
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FIG. 6: (Color online) The difference, A^xc, between the XC 
potential calculated at finite electric field and in absence of 
the field as a function of the site index. Results are presented 
for a 60 site chain with N = 10 (n = 1/6). The dots are the 
calculated data while the lines are a guide to the eye. The 
external potential has a negative slope. 
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FIG. 7: (Color online) The difference, A^xc, between the XC 
potential calculated at finite electric field and in absence of 
the field as a function of the site index. Results are presented 
for a 60 site chain with N = 30 (n = 1/2). The dots are the 
calculated data while the lines are a guide to the eye. The 
external potential has a negative slope. 



and reflecting the same trend already observed for the 
polarizabilities. 

There is however one anomaly in the potential response 
for the BALDA/LSOC functional, namely at n = 1/2 
and for small U/t (respectively 2 and 4) the potential 
response is actually opposite (positive slope) to that 
of the DMRG benchmark. This means that in these 
particular range of filling and interaction strength the 
BALDA/LSOC potential erroneously opposes to the ex- 
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ternal perturbation. The anomaly originates from the 
particular shape of the BALDA/LSOC potential as a 
function of n for small U/t (see figure [lj). In fact, v xc 
for BALDA/LSOC has a minimum for both U/t = 2 and 
U/t = 4 at around n = 1/4, which means that its slope 
changes sign when the occupation sweeps across n = 1/4. 
Therefore for those critical interaction strengths the re- 
sponse is expected to be along the same direction of the 
external potential for n < 1/4 and for 3/4 < n < 1 and 
opposite to it for 1/4 < n < 3/4 (at n ~ 3/4 there is a 
second change in slope). 

In the case of the BALDA/FN functional such an 
anomaly is in general not expected, except for small U/t 
and n close to the discontinuity at n = 1 (see figure Rl). 
This, however, is in the range of occupation not investi- 
gated here. Nevertheless we note that for n = 1/2 and 
U/t = 2 the BALDA/FN v X c is almost flat. This feature 
is promptly mirrored in the potential response of figure [7j 
which also shows an almost flat Awxc > although still with 
the correct negative slope. 

Given the good agreement for both the polarizability 
and the potential response between the exact results and 
those obtained with the BALDA (in particular with the 
FN flavour), one can conclude that the local approxima- 
tion to the Hubbard-LDFT functional is adequate. Still 
it is interesting to assess whether the remaining discrep- 
ancies have to do with the particular local parameteri- 
zation of Exc[{ n i}}, or with the fact that the exact XC 
functional may be intrinsically non-local. In order to an- 
swer to this question we have set a numerical test. We 
consider a 60 site chain with n = 1/2 (this should be long 
enough to resemble the infinite limit) and we introduce a 
local perturbation in half of the chain. This is in the form 
of a reduction of the on-site energy of the first 30 sites by 
S. We then calculate the deviation of the XC potential 
Sv as a function of the deviation of the total energy SEq. 
These two quantities are defined respectively as 

Sv = J2 l«xo(«<) - 4cV0l , SEo = E (S) - E o (0) , 

(16) 

with v xc and Eq (5) respectively the XC potential at site 
i and the total energy calculated for 6^0. One then 
expects for a local potential that Sv — > as 5E — » 0. 

Our results are presented in figure [8] These have been 
obtained for a relatively small U/t — 2 by varying S in the 
range < S < 0.1 in steps of 10 -5 (this range is used only 
for small <5, while a coarse mesh is employed for large S). 
Interestingly we note that, after a steady decrease of Sv 
with reducing SEq, the deviation of the potential starts 
to fluctuate independently on the size of SEq. We have 
carefully checked that such fluctuations are well within 
our numerical accuracy, so that they should be attributed 
to the breakdown of the local approximation. We then 
conclude that part of the failure of BALDA/FN in de- 
scribing the polarizability of finite ID chains must be 
ascribed to the violation of the local approximation. 




FIG. 8: Variation of the XC potential, Sv, as a function of the 
variation of the total energy, SEo, for a 60 site chain in which 
the first 30 sites have an on-site energy lower by 8 with respect 
to the remaining 30. The variation are calculated with respect 
to the homogeneous case. The inset shows a magnification of 
the data for small SEo. 



V. CONCLUSIONS 

In conclusion, we have reported a systematic study 
of the electrical response properties of one-dimensional 
metals described by the Hubbard model. This is solved 
within LDFT and local approximations of the exchange 
and correlation functional. Whenever possible the cal- 
culations are compared with exact results obtained ei- 
ther by exact diagonalization of with the density matrix 
renormalization group approach. In general we find that 
BALDA functionals perform rather well in describing the 
electrical polarizability of finite one-dimensional chains. 
The agreement with exact results is particularly good in 
the case of numerically evaluated functionals. A similar 
good agreement is found for the exchange and correlation 
potential response. In this case we obtain the interest- 
ing result that the potential response is always along the 
same direction of the perturbing potential, in contrast to 
what happens in ab initio DFT. Furthermore for small 
electron filling and weak Coulombic interaction the com- 
monly used LSOC parameterization is qualitatively in- 
correct due to a spurious minimum in the potential as a 
function of the site occupation. Finally we provide a nu- 
merical test of the breakdown of the local approximation 
being the source of the remaining errors. 
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